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Abstract. Applying the mean surface density of compan- 
ions, £(r), to the dynamical evolution of star clusters is 
an interesting approach to quantifying structural changes 
in a cluster. It has the advantage that the entire den- 
sity structure, ranging from the closest binary separations, 
over the core-halo structure through to the density dis- 
tribution in moving groups that originate from clusters, 
can be analysed coherently as one function of the stellar 
separations r. 

This contribution assesses the evolution of £(r) for 
clusters with different initial densities and binary popula- 
tions. The changes in the binary, cluster and halo branches 
as the clusters evolve are documented using direct A^-body 
calculations, and are correlated with the cluster core and 
half-mass radius. The location of breaks in the slope of 
£(r) and the possible occurrence of a binary gap can be 
used to infer dynamical cluster properties. 

Key words: binaries: general - open clusters and associ- 
ations: general - stars: formation - stellar dynamics 



1. Introduction 

Studying the clustering properties of stars in star-forming 
regions is a necessary input towards understanding the for- 
mation and evolution of young stellar clusters (Elmegreen 
& Efremov 1997). Using two-point angular correlation 
functions to analyse the spatial distribution of stars, 
Gomez et al. (1993) showed that young stars in the 
Taurus- Auriga molecular cloud form in small associations 
containing of the order of ten stellar systems. Larson 
(1995) extended this investigation by taking into account 
the results from different searches for binary companions 
to the pre-main sequence stars in the Taurus- Auriga re- 
gion. He computed the mean surface density of compan- 
ions, £(0), per star as a function of angular separation 
6. This statistical measure is closely related to the two- 



Send offprint requests to: R. Klessen 



point correlation function but does not require normalisa- 
tion. As two different power laws are necessary to fit the 
data, with a slope ks —2 for separations below rs 0.04 pc 
and with a slope ks —0.6 above, Larson concluded that 
there are two distinct clustering regimes. At small sepa- 
rations, the derived companion density is determined by 
binaries and higher-order multiple systems, whereas at 
large separations the overall spatial structure of the stel- 
lar cluster is observed. Larson interpreted the observed 
non-integer slope in the clustering regime as evidence for 
fractal structure of the cluster. In addition, Larson noted 
that the break of the distribution occurs in Taurus- Auriga 
at length scales which are equivalent to the typical Jeans 
length in molecular clouds, corresponding to a Jeans mass 
of w 1 M Q . He speculated that stellar systems with smaller 
separations form from the fragmentation of single collaps- 
ing proto-stellar cores, whereas the spatial distribution on 
larger scales is due to the hierarchical structure of the 
parent molecular cloud. 

These results prompted the reanalysis of the Taurus- 
Auriga data by other authors (Simon 1997, Bate, Clarke 
& McCaughrean 1998, Gladwin et al. 1999), as well as 
the subsequent investigation of additional star-forming re- 
gions. Mean surface densities of companions as a function 
of angular separation have been derived for Orion (Simon 
1997, Bate et al. 1998, Nakajima et al. 1998), for the p- 
Ophiuchus cloud (Simon 1997, Bate et al. 1998, Nakajima 
et al. 1998, Gladwin et al. 1999), and for the star- forming 
regions in Chamaeleon, Vela and Lupus (Nakajima et al. 
1998). Common to all studies is that the companion sur- 
face density is best described as a double power-law, with 
slopes of w —2 in the binary branch and slopes between 
—0.9 and —0.1 in the large-scale clustering regime. How- 
ever, the length scales where the break of the distribu- 
tion is found to vary considerably, from rs 400 AU in the 
Trapezium cluster in Orion, over k, 5000 AU in Ophiuchus 
and » 8000 AU in Taurus, to 30000 AU for the Orion OB 
region. This fact raises considerable doubts about the in- 
terpretation of the break location as being determined by 
the Jeans conditions in the cloud. This would imply quite 
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different Jeans masses which in turn should lead to devi- 
ations of the initial mass function, which have not been 
observed. 

A thorough theoretical evaluation of the mean sur- 
face density of companions, £(0), and a discussion of vi- 
able interpretations can be found in Bate et al. (1998). 
Altogether the following picture emerges: At small sepa- 
rations, £(0) traces the separations of binary stars and 
higher-order multiple stellar systems. The slope « —2 
results from the frequency distribution of binary separa- 
tions being roughly uniform in logarithm (Duquennoy & 
Mayor 1991, for main sequence stars). The break occurs 
at the 'crowding' limit, i.e. at separations where wide bi- 
naries blend into the 'background' density of the cluster. 
At larger separations, £(#) simply reflects the large-scale 
spatial structure of the stellar cluster. Bate et al. (1998) 
pointed out that £(#) can be strongly affected by bound- 
ary effects and that a non-integer power-law slope in the 
cluster branch does not necessarily imply fractal struc- 
ture. They showed that a simple core-halo structure, as 
is typical for evolved stellar clusters, will result in a non- 
integer slope for separations larger than the core radius. 
They also speculated about possible effects of dynamical 
cluster evolution on the properties of £(0). 

It is the aim of the present paper to investigate, for the 
first time, evolutionary effects on £(0) as derived from re- 
alistic TV-body computations. We use models studied by 
Kroupa (1995a, b,c, 1998, hereinafter K1-K4) for a com- 
parison with a 'standard' dynamical analysis, where the 
binary population is analysed separately from the bulk 
cluster properties. For comparison, the mean surface den- 
sity of companions for models of protostellar clusters that 
form and evolve through turbulent molecular cloud frag- 
mentation is discussed in Klessen & Burkert (2000, 2001). 

The structure of the paper is as follows. In the next 
section (§|J) we mathematically define the mean surface 
density of companions, £(r), and briefly discuss its limita- 
tions. In we describe the star cluster models and their 
properties. In §[|, we investigate the influence of cluster 
evolution on £(r) and in particular discuss wide-binary 
depletion. The effects of averaging and projection are anal- 
ysed in and a possible observational bias is discussed 
in §^|. §0discusses features in £(0) and their relation to 
cluster morphology. Finally, our results are summarised in 



2. Mean Surface Density of Companions 

The mean surface density of companions, £(0), specifies 
the average number of neighbours per square degree on the 
sky at an angular separation 9 for each cluster star. Know- 
ing the distance of the cluster, the angular separation 9 
between two stars in the cluster translates into an absolute 
distance r. In the numerical models we have full access to 
all phase space coordinates and we define, in what follows, 
the mean surface density of companions £(r) as the num- 



ber of stars per pc 2 as a function of the projected distance 
r (in pc). In Section 5.2 we show that the results are in- 



fer each star i in the system the projected distance 
all other stars j i. The separations 



variant to which plane is used for projection. We calculate 

to 

are sorted into 
annuli with radii r and width 5r, where we use logarithmic 
binning such that each decade in separation is partitioned 
into 10 logarithmically equidistant bins (log 10 Sr = 0.1). 
To obtain the function £(r), we divide the number SN(r) 
of stellar pairs per annulus r by the surface area 2irr5r 
and average by dividing by the total number N of stars in 
the cluster. The mean surface density of companions as a 
function of separation r then follows as 



E(r) = SN U 



2irNr5r 



(1) 



The function £(r) is related to the two-point correlation 
function, f (r), by f (r) = (£(r)/(£)) - 1, where (£) is the 
mean stellar surface density in the considered area (Pee- 
bles 1993). Because the normalisation (£) is often difficult 
to determine, it is preferable to use the function £(r). 

Stellar surveys have finite area and boundary effects 
may occur. For stars closer than a distance n, to the 
boundary all annuli with r > r& extend beyond the limits 
of the survey and companion stars may be missed. This 
has little effect when considering small separations r as 
only a small fraction of all stars is affected, however, when 
r becomes close to the survey size the missing companions 
result in a steep decline of £(r). Several methods have 
been proposed to correct for that effect (Bate et al. 1998, 
and references therein), where either the accepted range 
of r is reduced to separations much smaller than the sur- 
vey size, or additional assumptions about the background 
density of stars are made. Neither approach is completely 
satisfactory. 

In the present study we do not attempt to adopt any of 
the correction methods for the following reasons. First, in 
the numerical simulations all information about the sys- 
tem is accessible. There are no observational constraints 
and we always consider all stars in the cluster. The sur- 
vey area can be arbitrary large and is chosen such that it 
includes the complete cluster. Second, although the con- 
sidered clusters are subject to the tidal field of the Galaxy 
we do not include Galactic-field stars in consideration. 
Hence, there is no confusion limit, where contamination 
with foreground or background stars becomes important. 



3. Star Cluster Models 

A range of star-cluster models are constructed, and their 
dynamical evolution is calculated using Nbody5 (Aarseth 
1999), which includes a standard Galactic tidal field (Ter- 
levich 1987). The cluster models have been discussed ex- 
tensively in K1-K4, so that only a short outline is provided 
here. 
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Model A 


Model B 


Model C 


Model V 


Model £ 


Model T 


n 




5 


5 


5 


5 


3 


3 


f 




1 


1 


1 


1 










(pc) 


0.08 


0.25 


0.8 


2.5 


0.08 


0.25 


Th 


(10 6 yr) 


0.094 


0.54 


3.0 


17 


0.094 


0.54 


logio Po 


(stars/pc 3 ) 
(stars/pc 2 ) 


5.6 


4.1 


2.7 


1.1 


5.6 


4.1 


logio E 


4.3 


3.3 


2.3 


1.2 


4.3 


3.3 


t = tinit 


(10° yr) 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


t = 12.5 r h 


(10 6 yr) 


1.19 


6.67 


37.5 


211 


1.19 


6.68 


t = 125 r h 


(10 6 yr) 


11.9 


66.7 






11.9 


66.8 


t - ^end 


(10 6 yr) 


297 


300 


300 


296 


297 


300 



Table 1. Properties of the model clusters. All clusters have N = 400 stars with a mean stellar mass of 0.32 M Q , n 
lists the number of different realizations of the model, and / gives the starting binary fraction. Initially, the density 
profile of all clusters follows a Plummer law with half-mass radius rh, leading to a half-mass diameter crossing time 
Th, and to central and central surface densities po and So, respectively. The times, when we analyse the clusters, are 
indicated in the last four lines. 



The stellar systems initially follow a Plummer density 
distribution (Aarseth, Henon & Wielen 1974) with half- 
mass radius rh, and the average stellar mass is indepen- 
dent of the radial distance, r, from the cluster centre. 

Stellar masses are distributed according to the solar- 
neighbourhood IMF (Kroupa, Tout & Gilmore 1993) with 
0.1 < m < 1.1 Mq. Larger masses are omitted so as to 
avoid complications arising from stellar evolution. Binaries 
are created by pairing the stars randomly, giving a birth 
binary proportion / = N hin /(N sing + N hin ), where 7V sing 
and TVbin are the number of single-star and binary sys- 
tems, respectively. The initial mean system mass is 2(m), 
with (to) = 0.32 M Q being the average stellar mass. This 
results in an approximately flat mass-ratio distribution 
at birth (figure 12 in K2). Periods and eccentricities are 
distributed following Kl. The initial periods range from 
10 3 to 10 7 ' 5 days, and the eccentricity distribution is ther- 
mal, i.e. the relative number of binaries increases linearly 
with eccentricity being consistent with observational con- 
straints. 

The parameters are listed in Table |l|. Four clusters with 
/ = 1 are constructed spanning a wide range of central den- 
sities, from logioPo = 1-1 to 5.6 [stars/pc 3 ]. Each model 
contains N — 400 stars and has a mass M c \ = 128 M Q . 
The initial tidal radius in all cases is r t ~ 8 pc. All stars are 
kept in the calculation to facilitate binary-star analysis, 
but those with r r t experience unphysical accelerations 
in the linearised local tidal field and rotating coordinate 
system (Terlevich 1987), so that the density distribution 
of stellar systems at large radii does not reflect the true 
distribution in the moving group. Five different renditions 
are calculated for each model to increase the statistical sig- 
nificance of our results. In addition, two clusters with / = 
are constructed for comparison with the binary-rich cases. 
The evolution of these models is calculated for three clus- 
ter realisations each. The computations cover 1 x 10 9 years, 



but we consider only a sub-set of all possible snap-shots 
in the current analysis. 

4. Evolutionary Effects 

In this section, we discuss the influence of the dynamical 
cluster evolution on the resulting mean surface density of 
companions £(r). For each star cluster, £(r) is calculated 
as an average over the set of n individual model realisa- 
tions, and we restrict ourselves to discussing the projection 
into the xy-jA&ne. The influence of averaging and projec- 
tion is discussed in §0. 

4-.1. The Global Evolution of Stellar Clusters 

The evolutionary sequence of £(r) for models A to T> is 
illustrated in figure [l|. Its first column denotes the ini- 
tial state of each system at t = 0, the second and third 
columns depict the function £(r) taken at t = 12.5 Th 
and t = 125 Th, where Th is the initial half-mass diameter 
crossing time (table 1). As Th increases with increasing rh, 
these equivalent stages of cluster evolution correspond to 
different absolute times as indicated in the figure. The fi- 
nal state of the systems, at roughly t«3x 10 8 yr, is given 
in the last column. For models C and T>, 125 Th > £ C nd, so 
that the third column in figure [l] contains no entry. 

At any time and for all models, the mean surface den- 
sity of companions exhibits, at low separations r, a well 
defined power-law behaviour £(r) oc 1/r 2 . It implies an 
approximate uniform distribution of binary separations in 
logr (Bate et al. 1998). For larger r, the binary branch 
blends into the plateau of constant companion density 
corresponding to the core of the star cluster. This first 
break of the distribution occurs at separations ri where 
the number of chance projections of cluster members be- 
comes equal to the number of binaries in that separation 
bin. At larger r chance projections completely dominate 



4 Ralf Klessen and Pavel Kroupa: The mean surface density of companions 



t = 0.0 t=12.5r 1/2 t=125T 1/a t = t end 




Fig. 1. The mean surface density of companions, X(r), as a function of separation r for star clusters A to T> with 
an initial binary fraction of 100% at different states of the dynamical evolution: initially (t = 0.0, left column), at 
t = 12.5 Th and t = 125 Th (2. and 3. column) and at t = t cn( j (right column). The corresponding time in units of 
10 6 years is indicated in the upper right corner of each plot. S(r) is obtained as an average over n = 5 different cluster 
realizations for each model as a projection into the xy-plane. The error bars indicate Poisson errors. 



S(r), and wide binaries, if in fact present, can no longer be 
identified as such. The location of the first break therefore 
depends on the binary fraction and on the central density 
of the cluster. On larger scales the stellar clusters follow a 
Plummer radial density profile with half-mass radius rv 
Therefore, a second break occurs at r-i and S(r) declines 
sharply for separations r ^ w T2 (see Section 0). 

As the dynamical evolution progresses, the clusters ex- 
pand and the density declines. Hence, the projected mean 
surface density of companions decreases as well. While 
many of the binaries with separations comparable to the 
mean distance of stellar systems in the cluster core become 



disrupted, some new binaries may form by capture. Usu- 
ally these are higher-order multiples (K2) with separations 
close to the first break or smaller. As the cluster expands, 
the binary branch becomes less affected by crowding and 
the first break in E(r) shifts to greater separations. This 
behaviour is clearly visible in figure [l]. For all models the 
core plateau in E(r) 'decreases in height' and 'moves' to 
larger separations as time progresses. At late stages of the 
evolution the entire binary branch is uncovered and a few 
long-period orbits appear through capture. This is also 
documented in figs. 3 and 4 in K4. 
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300 
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Fig. 2. The mean surface density of companions, S(r), as a function of separation r for star clusters Z and J 7 . These 
models are equivalent to A and £>, respectively, but initially contain no binary stars. The figure is analogous to figure 
|l|, except that n — 3 different realizations of each model are used in the averaging process. 



The clusters develop core-halo structures through en- 
ergy equipartition. Low-mass stars gain kinetic energy 
through encounters with more massive stars. The low- 
mass stars move away from the cluster centre, forming 
the halo, whereas the more massive stars sink towards 
the centre. The trajectories of halo stars that trespass be- 
yond the tidal radius of the cluster are dominated by the 
Galactic tidal field, and most become unbound. Hence, 
the clusters expand until they fill their tidal radii. When 
this stage is reached (roughly after t > 1 x 10 8 years), 
the different clusters evolve identically (see also figure 1 
in K3), mostly loosing stars through their first and sec- 
ond Lagrange points (Terlevich 1987; Portegies Zwart et 
al. 2001). As a result, £(r) extends beyond the second 
break to increasingly larger separations with an increas- 
ingly shallower slope, which we identify as the halo branch 
in S(r). Near the tidal radius r t , a third break occurs, as 
stars with r > r t become unbound. The trajectories of 
stars belonging to these unbound moving groups are not 
followed with sufficient resolution (see Section |^). Also, 
these stars are likely to be severely contaminated by field 
stars in the Galaxy, and we refrain from a further discus- 
sion of this outermost branch in £(r). 

As can be seen from figure [l], knowledge of the ini- 
tial global properties of the system is effectively erased 
through the dynamical evolution. The cluster and halo 
branches in X(r) look quite indistinguishable in the final 
frames. The situation changes, however, w hen c onsidering 
the binary branch, as discussed in Section L2 . The / = 1 
versus / = experiments demonstrate that there is no sig- 



nificant difference in bulk cluster evolution between clus- 
ters containing a large primordial binary proportion and 
no binaries (K3). This is also evident by studying S(r); 
concentrating only on the cluster and halo branches, the 
upper two final panels in figures [j] and ^| are indistinguish- 
able. 



4-2. Binary Stars 

Binary systems and wide hierarchical systems form 
through capture during the evolution of the clusters, as 
seen in figure |^. These systems result from triple or higher- 
order stellar encounters, and are consequently very rare. 
The periods of these systems rang e from 10 7 to 10 11 days 
(figure 10 in K2), and a well distinguished binary branch 
with slope « — 2 develops, extending to radii well beyond 
the initial position of the first break. 

However, as the clusters evolve, binaries are not only 
created but also destroyed. Due to their smaller binding 
energies, wide binaries are more vulnerable to dynami- 
cal processes than close ones. If the initial binary fraction 
is high, then the destruction processes dominate over bi- 
nary formation, and cluster evolution leads to a depletion 
of wide binaries. As a result, S(r) steepens on the large- 
separation side of the binary branch. In extreme cases, 
some annuli r of S(r) may become completely depopu- 
lated, and consequently a gap between the binary and the 
cluster branches opens up, as is noticeable in figure |l|. 

The efficiency of binary disruption depends strongly 
on the initial density of the cluster. For high stellar den- 
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t = 0.0 t = 12.5 r 1/2 t=125r 1/2 t = t ei 




-6-5-4-3-2-1 -6-5-4-3-2-1 -6-5-4-3-2-1 

l°gio r lo gio r logio r 



Fig. 3. The function r 2 S(r) for models A to V. The depicted times, averaging and projection are analogous to figure 
[l]. The plots concentrate on the properties of the binary branch, where the horizontal dashed lines indicate its initial 
slope. The depletion of the distribution at late stages in the interval 10~ 4 pc < r < 0.1 pc is the result of wide-binary 
disruption. The rising part of r 2 at separations greater than the first break corresponds to the cluster core, and 
the following decline for r > n, is the contribution from the cluster halo. To demonstrate the effects of dynamical 
evolution, the dotted lines for t > 0.0 indicate the initial distribution. 



sities, the typical impact parameters of stellar encounters 
are small. Hence, there is a relatively high frequency of en- 
counters for which the energy exchange exceeds the bind- 
ing energy of typical binary systems, which subsequently 
dissolve. In our suite of models, the effect of binary deple- 
tion is largest in A, which has the highest central density 
Po, and decreases with increasing half-mass radius as 
Po becomes smaller. 

This is demonstrated in figure |^. Unlike the previous 
figures it shows a reduced range of separations, concen- 
trating on the binary branch, and it plots r 2 S(r) to make 



it easier to determine the power-law slope and deviations 
from it. At t = 0.0, for all models r 2 £(r) is constant in 
the binary branch. This reflects the initial distribution of 
the binary separations, which is uniform in the logarithm 
within the range 10~ 5 pc to 10~ 2 pc. Values r < 10 _5 pc 
come from the projection of the 3-dimensional distribu- 
tion into the xy-plane. For model T> the complete binary 
branch is visible and fully segregated from the cluster 
branch (which corresponds to the rising part of the plot). 
This is because the half-mass radius of the cluster is large 
enough that the projected mean separation between clus- 
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ter members exceeds the separation of the widest binary 
system. As both branches are clearly separated initially, 
dynamical evolution does not alter the binary distribution 
significantly. There is little sign of wide binary depletion, 
even at t = t cn d- 

In the other models, the typical separations in the 
cluster core are smaller than 10 -2 pc, and binary and 
cluster branches overlap in the beginning. This is most 
significant for model A, where the initial central den- 
sity is highest. Consequently a large number of wide bi- 
naries are disrupted during the dynamical evolution of 
the system, and r 2 E(r) drops considerably below its ini- 
tial value (indicated by the dashed line) in the range 
10 _4 pc< r < 10 _1 pc. Because the size of this gap de- 
pends on the age and the initial central concentration of 
the cluster, analysing the signatures in S(r) could be used 
to constrain the initial state of observed stellar clusters. 
This fact, namely that the binary population retains a 
memory of its past dynamical environment, is also used in 
Kl to infer the typical structures in which most Galactic- 
field stars form, by studying the shape of the binary period 
distribution ("inverse dynamical population synthesis"). 



5. Statistical Effects 

In the previous section we analysed the mean surface den- 
sity of companions averaged over several model realisa- 
tions seen in one projection (into the aiy-plane). Both, the 
effect of averaging and the effect of projection shall be 
discussed here. 

5.1. Cluster to Cluster variations 

To estimate the effect of statistical variations between dif- 
ferent model realisations, figure |J plots £i(r) constructed 
for each of the n = 5 individual cluster rendition of model 
A at t = 125 Th. This model is chosen, because it has 
the smallest initial crossing time, hence, our simulations 
span the largest evolutionary interval. After the system 
has evolved for 125 Th, differences become noticeable at 
separations where the number of neighbours is small. This 
is the case at very large separations, where the stellar den- 
sity rapidly decreases, and deviations may also occur at 
the extreme ends of the binary branch. At the smallest 
separations, S(r) is determined by only one or two very 
close binary systems, and near the first break of the dis- 
tribution, the depletion of wide binary systems becomes 
noticeable and a gap may open up. As can be seen in figure 
^, the exact location and number of depleted separation 
bins slightly varies between different model realizations. 
Therefore, the wide binary gap appears wider and more 
noticeable for individual realisations compared to the en- 
semble average. 



5.2. Different Projections 

Clusters of young stars are seen on the sky in only one 
projection. Inferring the full 3-dimensional structure of 
the cluster is therefore in principle impossible without ad- 
ditional information or assumptions. In the previous sec- 
tions, we concentrated on the projection into the xy-plane 
when plotting S(r). For the spherical clusters considered 
in the current analysis the different projections are equiv- 
alent, and each gives a fair representation of the complete 
system. Only at distances comparable to the tidal radius 
does the flattening through the Galactic tidal field of the 
cluster become significant (Terlevich 1987). 

The invariance to changes in projection is demon- 
strated in figure ||, which shows cluster .41 again at 
t = 125 Th. The function S(r) is essentially independent 
of the projection, and slight differences occur only where 
there are smaller numbers of stars. Analogue to the varia- 
tions between different model realisations discussed above, 
this is the case at very large and very small r, and at sepa- 
rations where wide-binary depletion occurs. At large sep- 
arations, the tidal field breaks the symmetry. For small r 
and at the binary gap, it depends on the projection which 
separation bin stays populated and which may become 
empty. However, besides these details the overall struc- 
ture of E(r) is projection invariant. 

6. Observational Bias 

The most likely observational bias in surveys of stellar 
clusters results from the unavoidable detection flux limit, 
which corresponds to a stellar mass limit m m i„ such that 
stars with m < m min are not detected simply because 
they are too faint. This effect influences the derived mean 
surface density of companions, as illustrated in figure |[ 
where we again concentrate on cluster model A. As in the 
previous figures, we plot S(r) at four different stages of the 
cluster evolution, however, we now 'observe' the cluster 
at various distances, i.e. we introduce different detection 
limits 7n m i n for the mass. We consider m m ; n in the range 
0.25 M Q to 1 Mq. 

With increasing m m i n , the total number of detected 
stars decreases, and as a result S(r) is reduced. Also the 
shape of E(r) changes. This effect is small for low detection 
thresholds (m ^ 0.5 M Q ), as the overall star distribution in 
the cluster is still well sampled. However, it becomes sig- 
nificant at large cluster distances when only the brightest 
stars can be detected. The binary branch becomes severely 
under-sampled, and wide gaps open up. For m m i n = 1 Mq 
the binary branch disappears in all models. 

The inferred cluster core radius also depends quite 
sensitively on the completeness of the stellar sample. As 
the clusters evolve dynamically, high-mass stars sink to- 
wards the cluster centre due to mass segregation, whereas 
low-mass stars move outwards, building up the extended 
halo (figure 2 in K3). Sub-populations of higher-mass stars 
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Fig. 4. Mean surface density of companions, S(r), of cluster A at time t = 125 Th for the projection into the xt/-plane. The 
figure illustrates the effect of the averaging process. The small plots show E(r) for the five different realizations Al to .45 of 
model A, whereas the large plot gives the resulting averaged function. 



therefore exhibit smaller core radii as time progresses rel- 
ative to the low-mass population. 

This trend is clearly seen in figure where the sec- 
ond break moves to smaller separations as m m i n increases. 
At very late stages of the dynamical evolution and for 
very large cut-off masses, the core radius may become 
too small, so that the second break is no longer notice- 
able. For example, the function S(r) follows an almost 
perfect r _2 -power-law for m m j n = 0.75 Mq at i C nd, ex- 
hibiting a smooth transition from the binary regime to the 
halo regime without any sign of the cluster core, which is 
present when taking all stars into account. When consid- 
ering only stars with m > 1 Mq , then the signature of the 
cluster core disappears at all times. This bias needs be 
taken into account when interpreting observational data 
on star clusters. 



7. Cluster Morphology 

As has been elucidated above, the mean surface density of 
companions shows distinct branches, the extend of which 
appear to couple with the dynamical state of the cluster. 
In this section we consider this in more detail. 



Simple bulk cluster properties that can be used to de- 
scribe the dynamical state of a cluster are the core radius, 
r c , the half-mass radius, rh, and the tidal radius, tv The 
core radius is approximated by calculating the density- 
weighted radius r c (Heggie & Aarseth 1992), 



^20 2 2 
1 'i Pi 



EN 2 o 2 
i=l l J \ 



(2) 



where pi = 3 771^5/(471" df 5 ) is the density around star i 
estimated within the closest distance di.5 from star i con- 
taining n — 5 additional stars with combined mass m;. 5, 
and r; is the distance of star i from the density centre 
of the cluster. The summation extends only over the in- 
nermost 20 % of all stars in the cluster, N20, since this 
is sufficient to ensure convergence. The approximate tidal 
radius (Binney & Tremaine 1987), 



nit) 



M st (t) 



3M 



gal 



(3) 



with Mgai = 5 x 10 Mq being approximately the Galac- 
tic mass enclosed within the distance of the Sun to the 
Galactic centre, tqc — 8.5 kpc. To estimate r t (t), M st (t) 
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Fig. 5. Illustration of the effect of projection. Mean surface density of companions E(r) for cluster realization .41 
at time t = 125 Th for three different projections. The upper panel shows the stellar distribution within ±2pc of the 
cluster centre and the lower panels show the resulting S(r). The global features of S(r) are independent of projection. 



is calculated by summing only those stars which have 
r(t) < 2r t (t — St op ), where the data output time interval 
St p <C ^reiax(i)- The quantities r c , and r% are averages 
of n models per time-snap (Table |l|). 

We define two breaks in E(r), n and ra, by fitting 
power laws to the three distinct branches of E(r) (the bi- 
nary branch, the flat central plateau, and the cluster halo 
out to the tidal radius r t ) and determining the separation 
at the intersection of the fits. These values are listed in 
Table ||, and plots of r c vs r% and vs are presented 
in figure [7| 

The figure shows that the quantities are well corre- 
lated. Specifically, we find (the uncertainties are mean ab- 
solute deviations) 

(r c /n) = 20.5 ±9.7 (4) 
and 

(r h /r 2 ) = 1.14 ±0.29. (5) 

This suggests that the break points in S(r) can be used to 
infer the core radius and half-mass radius of a cluster. We 
also find good correlations between and r c ((rh/r c ) = 
3.86± 1.04) and between the breakpoints in E(r) [r^/r-i = 
74.3 ±44.4). 

The above scaling relations apply for the specific low- 
mass cluster models that we investigate. These models do 
not suffer significant core collapse, which is partly given 



by the relatively fast evaporation time (« 0.5 — 1 Gyr, 
K3), and the ubiquitous binary stars which oppose core 
collapse. More massive clusters are likely to show different 
correlations, notably between and r c , and between r\ 
and rii since core collapse leads to the contraction of r c 
but an expansion of (Giersz & Spurzem 2000). Analysis 
of more massive clusters using E(r) is a future goal, and 
it will be interesting to see if correlations ^ and || remain 
valid. 

8. Summary 

In this paper, we investigate how dynamical cluster evo- 
lution is manifest in the mean surface density of compan- 
ions, E(r), as a function of separation r. We find that 
throughout all evolutionary phases, £(r) can be subdi- 
vided into four distinct branches, each following approxi- 
mately a power-law behaviour: 

1. At small separations, and throughout all evolutionary 
phases of the star cluster, E(r) traces binary stars and 
higher-order multiple systems. In the binary branch 
the slope of the companion density is approximately 
-2. 

The binary population keeps a memory of the dynam- 
ical evolution of star clusters. The depletion of the bi- 
nary population at large separations leads to a steep- 
ening of E(r) in the wide-binary branch. In extreme 
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Fig. 6. Mean surface density of companions, S(r), of cluster A for different observational minimum mass limits m m ; n . 
First row of panels: all stars in the cluster are considered (analogue to top row in figure ^). Second row: only stars 
with m > 0.25 Mq contribute to E(r). Third row: m m in = 0.50 M©. Fourth row: m m i n = 0.75 Mq, and lowest row: 
Wmin — 1.00 Mq. Times and averaging procedure are equivalent to figures I to | There are 400 stars in each cluster 
realization in the mass range 0.1 M© to 1.1 M© following the Kroupa et al. (1993) IMF. The fraction of cluster stars 
considered in each panel is 100%, 51%, 17%, 6%, and 1%, respectively. 



cases, it may become completely depopulated and a 
gap between binaries and the cluster branch opens up. 
The gap increases with cluster age and with increasing 
initial density. For the binary branch, E(r) is sensitive 



to the initial properties of the system, so that analysing 
S(r) for a cluster can constrain its birth configuration. 
2. The well-defined binary branch blends into a plateau of 
constant companion density which corresponds to the 
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Model T 
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0.25 
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Table 2. Morphologically important length scales. The table lists the separations r\ and where the first and second 
break of the distribution X(r) occur in all model clusters. The values for n and r 2 are obtained at the intersection 
of the power-law fits to the binary branch and the central flat plateau, and the plateau and the halo distribution 
of cluster stars out to the tidal radius, respectively. These separations are compared to the core radius r c and the 
half-mass radius rh at the different times in Fig. m. 



main body of the stellar cluster, and which we refer 
to as the cluster branch in S(r). The transition oc- 
curs at separations when chance projections of cluster 
members begin to outnumber the contribution of bi- 
nary stars in that separation bin. The location of this 
first break at r± depends on the binary fraction and 
the central stellar density of the cluster, and thus on 
r c . As the cluster expands, the density decreases and 
the core radius increases with r c w 21 r\ for the set 
of models studied here. As a consequence, the binary 
branch becomes less affected by crowding and the first 
break in S(r) shifts to greater separations. Generally, 
as the cluster evolves the plateau 'decreases in height' 
and 'moves' to larger separations. Larger parts of the 
binary branch are 'uncovered' as a result of the expan- 
sion and the formation of wide hierarchical systems 
through capture. 
3. Relaxed Galactic clusters exhibit a core/halo struc- 
ture. This becomes apparent in E(r) through the exis- 
tence of a second break at r 2 , which is approximately 
located at the half-mass radius, rh ~ r 2 , beyond which 
E(r) decreases again. The slope of this halo branch 
depends sensitively on the evolutionary state of the 
cluster. If the cluster is still sufficiently young so that 
most of it is confined well inside it's tidal radius, then 
E(r) decreases rapidly. However, relaxed clusters that 



fill their tidal radii have a halo branch with a slope of 
approximately —2. 

4. Beyond the tidal radius, the density of unbound stars 
decreases gradually, which is reflected in a fourth 
branch in E(r). It is a measure of the density distri- 
bution of the moving group relative to the centre of 
its origin, which is the star cluster. The decay of £(r) 
is steeper here because the stars are finally removed 
from the vicinity of the cluster within an orbital pe- 
riod about the Galaxy (Terlevich 1987). 

5. The mean surface density of companions closely re- 
flects morphological properties of stellar clusters. The 
break points in E(r), r\ and r 2 , therefore can be used 
to infer the core radius r c and half-mass radius n, of 
a cluster. These quantities are well correlated, and we 
find (r c /ri) = 20.5 ± 9.7 and (r h /r 2 ) = 1.14 ± 0.29, 
respectively. 

Altogether, £(r) contains valuable information on the 
dynamical state of a star cluster, and allows an assessment 
of cluster properties such as the core and half-mass radii. 
As the location of the first break in £(r) (or the possible 
occurrence of a gap between binary and cluster branch) 
depends on the binary fraction and on the initial central 
density, these parameters can in principle be inferred from 
analysing E(r). 
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Fig. 7. Correlation between the first break radius n and cluster core radius r c (left panel), and the radius ri of the second 
break and the half mass radius r h (right panel). 



Our study confirms that different projections of the 
same data do not change S(r) significantly during the 
evolution of initially spherical clusters in the Galactic tidal 
field. Also, different numerical renditions of the same mod- 
els lead to indistinguishable results. Hence, they can be 
combined to improve the statistical significance of the en- 
semble average S(r). 

To allow for a proper assessment of stellar cluster prop- 
erties using £(r), it is important to consider as complete a 
census of cluster stars as possible. Observational bias (i.e. 
non-detection of faint stars) may complicate the interpre- 
tation of £(r) and limit its applicability for inferring the 
initial stage of the cluster under study (see also Bate et 
al. 1998). Our tests show, however, that E(r) remains a 
useful quantity even when stars with m^0.5M a are not 
detected. We also find that mass segregation is evident in 
E(r) through the location of the second break in depen- 
dence of the mass-range used to construct S(r). 

Future analysis of numerical models of rich clusters, 
for which mass segregation and possibly core collapse play 
important roles in the late phases of the dynamical evolu- 
tion, will be performed to deepen the issues raised in this 
pilot study. 
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